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Abstract. We investigate the unification scenario provided by the generalised Chaplygin gas 
model (a perfect fluid characterized by an equation of state p = —A/p°'). Our concerns lie 
with a possible tension existing between background kinematic tests and those related to the 
evolution of small perturbations. We analyse data from the observation of the differential age 
of the universe, type la supernovae, baryon acoustic oscillations and the position of the first 
peak of the angular spectrum of the cosmic background radiation. We show that these tests 
favour negative values of the parameter a: we find a = -0.089l:[]:|^^ at the 2a level. These 
would correspond to negative values of the square speed of sound which are unacceptable 
from the point of view of structure formation. We discuss a possible solution to this problem, 
when the generalised Chaplygin gas is framed in the modified theory of gravity proposed by 
Rastall. We show that a fluid description within this theory does not serve the purpose, 
but it is necessary to frame the generalised Chaplygin gas in a scalar field theory. Finally, 
we address the standard general relativistic unification picture provided by the generalised 
Chaplygin gas in the case a — )• 0: this is usually considered to be undistinguishable from 
the standard ACDM model, but we show that the evolution of small perturbations, governed 
by the Meszaros equation, is indeed different and the formation of sub-horizon GCG matter 
halos seems to be strongly suppressed in comparison with the ACDM scenario. 
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1 Introduction 

Under the assumption of homogeneity and isotropy of the universe at large scales, and that 
General Relativity is the correct theory of gravity, the accelerated expansion [1-4] is due to 
an exotic component, with negative pressure, called dark energy (DE). Moreover, structures 
such as galaxies and clusters of galaxies seem to require the existence of another exotic 
component, dubbed dark matter (DM) [5, 6], in order to form. DE and DM constitute the 
so-called dark sector of the universe, accounting for about 95% of the total cosmic energy 
budget. Their nature is still a mystery, but they clearly do not fit in the standard model of 
elementary particles [5, 6]. 

One important proposal describing the dark sector of the universe is the unification 
scenario, where DM and DE are considered diff'erent manifestations of a unique component. 
The paradigm of such idea is the Chaplygin gas (GC) model [7], where an equation of state 
inspired by string theory [8] allows to obtain a dynamic behaviour mimicking DM in the 
past and DE in current time, as structure formation and the accelerated expansion require. 
Later, the CG has been phenomenologically generalized, leading to the so-called generalized 
Chaplygin gas (GCG) [9, 10]. 

Though the idea is very appealing, the GCG faces crucial drawbacks when the model 
undergoes observational tests. This has been first indicated in reference [11], concerning the 
behaviour of small perturbations and structure formation, and pointed out again in [12-17]. 
One important aspect is that the so-called kinematic tests (those based on the background 
expansion and the calculation of distances) seem to favour values for the GCG equation of 
state parameters that imply negative square speed of sound [18, 19] (though they do not 
exclude positive values). If this is really the case, up to our knowledge in the literature there 
are only two ways out: i) to introduce ad hoc entropic perturbations [20, 21] which nullify (or 
make positive) the effective speed of sound or ii) to modify the gravity theory, by supposing 
e.g. a non-conservative theory of gravity [22]. 
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Our goal in this work is to understand to which extent the results of [18, 19] are really 
robust from the point of view of the kinematic tests. This question has already been addressed 
in [19], using the differential age of the universe (i.e. the H{z) test) and type la Supernova 
(SNIa). We extend that analysis by enlarging the number of observational data for the 
H(z) test and by including the position of the first acoustic peak of the cosmic microwave 
background (CMB) spectrum and the baryonic acoustic oscillation (BAO) tests. We revisit 
the constraints placed by these datasets, but we do not touch delicate issues such as the 
calibration problem of the SNIa data [23] . Our results seem to reveal the robustness of the 
kinematic tests for the GCG model: we find a = -Omrtoill 

(at the 2a level). Therefore, 

we not only confirm the previous results, but stress that negative values of a seem to be 
favoured over the positive ones (for comparison, see Table 1 of [19]). 

Note that we cannot report a real "tension" between background tests and the ones 
involving perturbations, since at 2a we find that positive values of a are still allowed (a < 
0.072, whereas in [11] the authors find |a| < 10"^, for perturbations), but we can speculate 
how the Chaplygin gas could be "saved" if such tension would appear in the future, given 
the fact that many experiments are ongoing and others forthcoming (such as EUCLID), thus 
the precision on the observational data is increasing day after day. 

Within the "tension" scenario, if the GCG is indeed the matter component behind 
the dark sector, it must get its origin from a different theory of gravity rather than GR. 
For example, non-standard gravity theories like Rastall [24-28] or f{R) [29] could possibly 
improve the status of the GCG. Our discussion in Sec. V is devoted to open a new window 
for GCG based cosmological models where standard GR is replaced by Rastall's theory. 

We also stress that, in general, the accepted idea is that the GCG (based on GR theory) 
works only if its parameter space is reduced to the ACDM one, in particular if a = 0, see 
Eq. (2.1) and Eq. (2.4) [30]. We discuss this issue in Sec. VI. We point out that the reduction 
a = implies that the GCG does not behave as a pressureless fluid at high redshifts, but as a 
fluid with a small negative constant pressure. Hence the analysis of the growth of sub-horizon 
GCG dark matter halos during the matter dominated epoch and the "Meszaros effect" could 
be different from the standard ACDM picture leading to a different clustering pattern at 
small scales and to a modification of the transfer function. We obtain a "Meszaros-like" 
equation for the GCG with a = and show that it does not reproduce the typical behaviour 
expected in the standard ACDM scenario. This argument agrees with the conclusions of [31]. 
For a opposite point of view, i.e. that the perturbations in the GCG with a = are the same 
as in the ACDM, see [32] . 

The paper is organized as follows. In Sec. 2, we review briefly the possible tensions 
within GCG model tests. In Sec. 3 we describe the observational tests we perform, while 
in Sec. 4 we carry out the statistical analysis based on such tests. In Sec. 5 we investigate 
a Rastall's theory approach which could possibly save the GCG. In Sec. 6 we discuss the 
general relativistic limit a — )• of the GCG and solve its Meszaros equation. In Sec. 7 we 
present our conclusions. 

2 Successes and problems of the Chaplygin gas model 

The GCG model is characterized by the equation of state [7, 9, 10] 

A , , 

^'c = --^, (2.1) 

rc 
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where A and a are free parameters. The original Chaplygin gas model, which is somehow 
connected with the Nambu-Goto action of string theory, implies a = 1 [8], the case a ^ 1 
being a phenomenological generalization. Integrating the conservation equation for the fluid, 

dpc 3 Pc -A 

— 1 — pc(l + Wc) = , with EoS parameter Wc= — = Yir~ i (2-2) 

da a Pc pc'^°' 

leads to the following expression for the GCG density as function of the scale factor: 

(1-^)- 



Pc{a) = PcO 



A + 



a- 



3(l+a) 



1 

l + Q 



(2.3) 



where we have defined the new parameter A = A/ p^'^^ and present-time quantities are 
indicated by the subscript 0. Prom Eq. (2.3) one can see that the original motivation behind 
the GCG, i.e. a fluid that evolves from the matter behaviour at early times to a constant 
density at late times, occurs only if a > — 1. 

When one consider a model where the energy content is given by radiation and baryons, 
besides the Chaplygin gas, Friedmann's equation becomes 



Hi 



Om(1 + zf + a,o(l + zf + OeO + (1 - ^)(1 + zf^^+^~^ 



1 

l + a 

, (2.4) 



OcO — 1 — ^bO — ^rO , (2.5) 

where Hq is the Hubble constant and z is the redshift, which is related to the scale factor 
by a = 1/(1 + z). We have also assumed a spatially flat background, according to the recent 
results of WMAP7 [3]. The Hubble parameter today can be expressed as Hq = lOO/i km 
s~^ Mpc~^, where h ~ 0.7 (we will adopt the range 0</i<<lasa prior for our Bayesian 
analysis). We also introduce E{z) = H{z)/Hq, for future convenience. 

The idea of DM-DE unification into a single fluid in expression (2.4), faces difficulties 
from different points of view. From the theoretical side, if vacuum energy is not responsible 
for DE, it should be explained how it contributes to gravity and therefore how it enters 
Eq. (2.4). From the observational point of view, introducing from Eq. (2.1) the GCG speed 
of sound 

A 

c^(a) = —awc{a) , if a = 1 — )■ c^q = a — tt = OiA , (2-6) 

PcO 

one can see that if background tests favour negative values of the parameter a, then 
is negative. This is not only an undesirable feature of the model, but depending on how 
much it is negative a "tension" between background and perturbative tests may arise. These 
problems can be alleviated if, for example, positive entropic perturbations a = dp/dS > are 
introduced such that the effective speed of sound becomes a positive quantity. However, the 
introduction of entropic perturbations implies new free parameters: there are some degrees 
of arbitrariness, and in some cases, such as those of reference [20], it is possible only a slight 
extension to negative a. Another possibility of achieving negative values of a is to implement 
a scalar version of the GCG model using a non-conservative theory of gravity, like Rastall's 
theory. This has been done in reference [22], and in fact negative values of a seem still to be 
favoured at the level of perturbations. 

This problem concerns not only the confrontation with the power spectrum data, but 
also the analysis of the anisotropy of the cosmic microwave background radiation (CMB). In 
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references [12, 13], the confrontation of the GCG model with the CMB spectrum has been 
performed, indicating that the most favoured scenario implies a — )• 0, that is, the GCG model 
reduces to the ACDM model. This result is one of the reasons for speculating that the GCG 
model should be be ruled out. 

3 The background tests 

There are four main background tests for a cosmological model: 

1. The differential age of old galaxies, given by H{z). 

2. The SN la data. 

3. The position of the first CMB acoustic peak. 

4. The peak position of the baryonic acoustic oscillations (BAO). 

For the differential age data, connected with the evaluation of the age of old galaxies 
that have evolved passively, leading to values of H{z) for specific redshifts, there are 13 
observational data [33-^37]. Recently, a compilation of 21 data points has been considered 
[38]. The fundamental relation is 

, ^ 1 dz , ^ 

The value of the Hubble parameter today obtained by the HST [Hq = 72 km s~^ Mpc~^) 
could also be added to this sample, but this would imply a prior on the final parameter 
estimation and, therefore, will not be included here. The analysis with the sample of 13 data 
points for the Chaplygin gas model has been made in reference [19], leading to the conclusion 
that slightly negative values for a are favoured. In our analysis below we use the set of 21 
data points compiled in ref . [38] . 

The SN la test is based on the luminosity distance, given by 



For the SN la data, we have two main problems. The first one concerns the choice of the 
sample. There are many different SN la data set, obtained with different techniques. In 
some cases, these different samples may give very different results. The second point is the 
existence of two different calibration methods: one using cosmology and which takes into 
account SN with high z (Salt2); the other using astrophysics methods, valid for small z 
(MLCS2k2) [23] . In some case, the employment of different calibrations can lead to different 
results also. All this, make the SN la analysis very delicate. Here, however, we use the 
Union sample [39], calibrated by the Salt2 method. This choice is motivated by looking for 
a contact with previous results in the literature, including those of reference [19]. 

The position of the first peak of the CMB spectrum is a more complex test. It is linked 
to oscillations of the baryon-photon plasma at the recombination period and is given by 



H = m - M = 51ogio-DL 



(3.2) 




(3.3) 



h=lA{l-(l)l) . 



(3.4) 
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A detailed numerical analysis leading to the following fitting formula for parameters of this 
fundamental quantity [40-43]: 

'Ai = 0-267 ( Y'\ (3.5) 

51 = 0.0783 (l + 39.5a;,V'') ' , (3-6) 
c,2 = 0.560 (^1 + 21. iLc-fe^o^^^ , (3.7) 
Zls = 1048 (^1 + 0.00124a;^-°-^38^ ^ ^^^92^^ ^ (3 8) 

= vr£ , (3.9) 

c-(3.^^(l..))""\ (3.11) 

We will use the following values for the different density parameter: 

ujio = nioh^, 0,.o/i^ = 4.116 X 10~5 , Ofeo/i^ = 0.02258 , comoh^ = ^dmo + ^bo ■ (3.12) 

From observation: h = 220.08 ± 0.7. 

The baryonic acoustic oscillations are due to the effect of oscillations in the photon- 
baryon plasma at the moment of the decoupling, at about z = 1090 (it is of course the 
same physics which produces the acoustic peaks structure in the CMB, but now this effect is 
observed in the baryonic distribution). This effect is quantified by the following expression 
[45]: 

vn;^ / I 



A = ^i- / ^ . (3.13) 

In this work we use data from the WiggleZ Dark Energy Survey [46]. 

The GCG behaves as dust in the past. Hence, we should identify an effective DM 
component in the GCG in order to use the above formulas for the analysis. There are different 
prescriptions in this sense in the literature. We will adopt the decomposition proposed in 
[44], where the effective DM component is given by, 

^mo = ^bo + nco(l-A) . (3.14) 



4 Analysis and results 

We perform a Bayesian statistic analysis in comparing the theoretical predictions with the 
observational data. First, for each dataset, we compute 

(4.1) 

1=1 
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where fif^ is the theoretical prediction for the quantity whereas is the correspond- 
ing observational datum, with an error estimation Uj. From the chi-squared, we build the 
probability distribution function (PDF), as follows: 

P{h,a,A) = Ce^p(^-^y (4.2) 

where we assumed the data to be independent and normally distributed; C is a normalization 
constant. As indicated, the PDF depends on three free parameters: h, a and A, characterizing 
a three dimensional function. Two and one dimensional PDF can be computed integrating 
(marginalizing) on the remaining parameters. 

As already said, we will compute the separately for each dataset. The total for 
the set of the four observational tests, shall be the sum of the separated x^- Accordingly to 
the exponential form of the PDF, in Eq. (4.2), the total PDF shall be the product of the 
single ones. 

The intervals considered for each of the free parameters are crucial for the final esti- 
mations. For A, all values < j4 < 1 are considered. For h we consider both delta and 
constant priors i.e., for the former prior we fix values for h while for the latter we integrate 
in the interval < h < 1. For a, the situation is more complex, since we can consider two 
possibilities: a > — 1, in order to assure a transition from dust in the past to a cosmological 
constant in the future, or a > 0, in order to assure a positive square speed of sound. Our 
goal here is to show that negative values for a are preferred. Therefore, we focus our analysis 
on the prior a > —1. However, for the sake of completeness, we also carry out a statistical 
analysis leaving a free to vary to arbitrarily negative and positive values. 

We evaluate the one dimensional PDF for A and a for the four independent tests listed 
before and combining all of them. The results are shown in Figure 1-6. 

4.1 Delta prior over h 

In Figure 1, 2 we show the one-dimensional probability distributions for a and A for specific 
choices of /i = . The values chosen for h* are shown in the four panels. Each panel 
corresponds to a data set. All the curves have been normalized following 



P{h\a,A)dA P{h\a,A)da 

—r^ — ^ _ and PDFM = ^ ' _^ 

£ f°° P{h\a,A)dadA £ f°° P{h*,a,A)dadA 



Except for the H(z) data, we see that the maximum likelihoods occur for hyper-surfaces with 
negative values of a. 



4.2 Flat prior over h 

For this analysis we left h to vary in the range < h < 1. The final one-dimensional PDF 
for a and A are calculated formally as 

fn P(h,a,A)dhdA nP{h,a,A)dhda 
PDF(a) = ^ _ and PDF(yl) = ^ . 

lo f~i /o ^(^' ^) '^^ lo /-I /o ^(^' ^) '^'^ 

(4.4) 
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It is crucial to restrict the interval for a. If we consider the standard scenario for structure 
formation, in the GR context without entropic perturbations, it is required that a > 0. 
However, if the only constraint is to impose an acceleration to deceleration transition, then 
the restriction is a > —1, since for a < —1 the universe accelerates in the past, and it 
is decelerating today. Alternatively, we may leave the value of a free in order to test the 
consistency of the background tests in the context of the GCG model. In what follows, we 
will consider two cases: a > —1 and a free. 

For both the cases, the general behaviour follows similar features, that can be summa- 
rized as follows. The H(z) and the baryonic acoustic oscillations tests predict a maximum for 
a slightly negative, while the position of the first peak indicates a maximum for a small pos- 
itive value of a. For these three tests the PDF decreases as the value a = — 1 is approached. 
Only for the SNIa test the PDF for a can be significant for a < —1. It is important to stress 
that, for the baryonic acoustic oscillations and for the CMB first peak, the PDF becomes es- 
sentially zero for a < —1. This is due to the decomposition into a "dark matter" component 
employed in Eq. (3.14): for a < —1, the behaviour for "dark matter" is reversed, due to the 
change of sign of the exponent in Eq. (3.14). 

We show the results of the analysis in Figure 3 and Figure 4, for a > — 1 , and in Figure 5 
and Figure 6, for a free. The 2a estimation for a in the first case is a = — 0.0871q'J|5, while 
in the second case is a = — O.OSOIq 'JIq. This result shows that a > —1 is a good prior. 



5 A Chaplygin gas fluid description in Rastall's theory and some possible 
indications for its healing 

Some ideas on how to avoid a possible "tension" and then saving the Chaplygin gas appeared 
recently and concern a modified theory of gravity: Rastall's theory [24]. In the latter the 
Einstein-Hilbert action still holds and the modification resides in the conservation law of the 
matter stress-energy momentum, i.e. 



where k is a parameter and R is the Ricci scalar curvature. The above modification can be 
interpreted in various ways. To us the most interesting and significant one is that Rastall's 
idea may be viewed as a kind of semi-classical formulation of quantum phenomena, which we 
expect to appear when the curvature (which enters as R on the right hand side of Eq. (5.1)) 
becomes important. Of course, one could choose other scalars rather than R in order to 
represent the Riemann tensor, but perhaps R is the most natural choice. 

Since Eq. (5.1) must fit into Bianchi identities, one can find the following modified 
Einstein equations 



T^^;^ = kR' 



(5.1) 




(5.2) 




7-1 



(5.3) 



2 



7 



Then, it is clear that, when assuming a FLRW background, Friedmann equation shall be 
modified with terms proportional to 7. 

We choose, for our discussion in the present section, a single-fluid component with 
density p and pressure p. In Rastall's theory, Friedmann equation becomes 



7 3i(; , ^, 
+ ^(7-1) 



2 2 



(5.4) 



where note that, differently from GR, now the pressure contributes to the expansion (via the 
equation of state parameter w). Moreover, the continuity equation becomes: 

p + 3F(p + p) = ^(p-3p) , (5.5) 

or, trading the cosmic time for the scale factor, and rearranging the derivatives: 

'i-l dp 3 2-1 dp 

Clearly, since the energy conservation equation has changed, we cannot use again the GCG 
equation of state and hope to obtain again the same result as in Eq. (2.3). What we can do 
is to take advantage of the foregoing analysis and assume a form for H given by Eq. (2.3). 
Then, we can derive the equation of state and speed of sound of the corresponding fluid in 
Rastall's theory and investigate its stability properties. 

Using E = H'^ /Hq, it is easy, by simple inspection of Eqs. (5.4) and (5.6) to find 

--{P + P) = ^, (5.7) 
a da 

which is a relation which holds true also for GR. Solving Eq. (5.4) for p and substituting in 
the above equation, one finds 

j-1 dE 1 , , 

-a-r + o — (5-8) 



and for the pressure 



Now, as said, we assume 



2(3 - 27) da 3-27 



'—E. (5.9) 



E 



6(3 - 27) da 3-27 
(1-^) 



A + 



7,3(1+") 



1 

\ + a 

(5.10) 



and then it is straightforward to compute density and pressure. In particular, the equation 
of state parameter and the speed of sound take the form 

,,^-jM;^±h^M^, (5.11) 

A (2a3"+3 + 3^ _ 5) _ 3^ + 5 ' ' 
2_ ^{7[(a + l)a3"+3-l] -(3Q + l)a3"+3 + l}+7-l 
'^^ ~ A {37 [{a + l)a3°+3 - 1] - (3q + 5)a3"+3 + 5} + 3^ _ 5 ' ^ ' ' 

where we have defined the speed of sound as = dp /dp = {dp /da) /{dp /da), i.e. as if it 
were adiabatic. We can formally invert the function p{a) and then obtain a p = p{p), i.e. a 
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barotropic equation of state. We plot in Figure 7 the above w and for the choice a = —0.1 
and A = 0.7 (the best fit values found in this paper) and for 7 varying about unity. 

In the range chosen for 7, the speed of sound is monotonically decreasing and, by 
inspection, positive if 7 ^ 0.8. Let us investigate briefly the asymptotic: in the far past, for 
a — )• 0, we have 



2 ,7-1 



37 - 5 ' 



, (5.14) 



and for 7 = 1 we reproduce the known GR result for the GCG. Now, in the remote future, 
for a — >■ 00: 

w ^ —1 , a — )• 00 , (5.15) 
2 , 7 - 1 + a(7 - 3) 



37 - 5 + 3a(7 - 1) 



00 , (5.16) 



where again we reproduce correctly the limit — >• a when 7 = 1. By construction, the final 
stage of the evolution is a de Sitter one. Of course, in all the above calculations we have 
assumed that a > —1. Asking for the asymptotic speed of sound of Eq. (5.16) to be positive, 
we find the following ranges for 7: 

3a + 1 5 — 3a , „s 

7 < -r-— , 7 > r • 5.17 

1 + a 3(1 — a) 

For a = — 1, this amounts to 7 < 7/9 and 7 > 53/33. 

5.1 Metric perturbations and evolution of small perturbations 

We now check the evolution of perturbations. Let us start by considering the energy- 
momentum tensor 

Tl = pu^u^ + , (5.18) 

where /i^*^ = g^^ + u^^u^ . More explicitly, the background components of (5.18) are 

A 

J- 

We assume a conformal Newtonian gauge line element for scalar perturbations 

ds^ = c? (77) [- (1 + 2$) (ir?2 + (1 - 2^) bijdx'dx^ , (5.20) 

where we introduced the conformal time r]. The perturbations of the fluid 4-velocity up to 
first order are given by 

«° = ^(l-$), 'Uo = -a(l + $). (5.21) 

In the absence of anisotropic stresses the spatial off-diagonal Einstein equation implies $ = 
Assuming the above perturbed metric, we can find the Einstein equation for the potential <^ 



T\ = -p, Tj^ = T^ = 0, T;=p6i = -^6). (5.19) 
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for the Rastall's theory, namely 

A$ - 3 



a! (a! 



$" + 3— ^>' + 

a 



a \a 



, a \ I a 
2(- + - 

a J \ a 



Sp — — {6p — 3Sp) 



7-1 

Sp H — {6p + 6p — 36p) 



(5.22) 



(5.23) 



where the prime denotes derivation wrt the conformal time. Introducing the speed of sound 
(assuming adiabaticity, i.e. 6p/6p = dp /dp): 



A^>-3— (—^ + ^' 



^>" + 3— ^>' + 

a 



a \ a 
a J \ a 



4^Ga^Spi^ + 3^c^ 



(5.24) 
(5.25) 



Changing the conformal time for the scale factor and combining the two above equations, we 
find an equation for the gravitational potential: 



n 



1 



Ua 



7 - 1 + (5 - 37)c^ 1 
3 - 7 + 3(7 - l)c2 a? 



$ 



3$ - 3a$a 



(5.26) 

where the subscript a denotes derivation wrt the scale factor, % = a' /a and we have intro- 
duced a plane- wave expansion. Note the "Rastall factor" of Eq. (5.26): 

1 + (5 - 37)c2 



Rf 



7 ■ 



Del 



(5.27) 



-7 + 3(7 

Since it multiplies the wavenumber k, it may be considered as an "effective" speed of sound, 
different from the adiabatic we introduced. In this sense, Rastall's theory seems to call 
into play a sort of "geometric entropy". Using Eq. (5.12) in order to reduce the "Rastall 
factor" one finds 

Aa 

= ^ + (i_^)a-3(a+i) ' (5-28) 

which is exactly the GCG speed of sound in GR! Rastall's theory seems to be able to reproduce 
the same evolution of perturbations of a fluid in GR which provides a given background 
expansion. Therefore, we may conclude that even framing the GCG in Rastall's theory does 
not save the model from being ruled out due to the behaviour of small perturbations. 

Just to check. We choose as initial conditions $ = — 1 and <I>a = in a = 0.005 and let 
the system evolve with the choice a = —0.1 and A = 0.7. In Figure 8 we plot the evolution 
of $ and 6 as function of the scale factor for a representative scale /c = 0.1 h Mpc~^. 

Clearly, the results are catastrophic. A remaining possibility of salvation is to adopt a 
scalar field description for the GCG, as done in [22]. As shown in [50] (where the authors do 
not mention Rastall's theory and investigate just the 7 = 2 case) and [27] in a scalar field 
approach the effective speed of sound in the rest frame of the field (i.e. where Tq = 0) can 
be written as 

2_ 2-7 



7 



(5.29) 



which evidently vanishes for 7 = 2. Since the scalar field description possesses one degree of 
freedom more than the fluid description (encoded in the scalar field potential) the background 
expansion can be fixed without any drawbacks on the evolution of perturbations, which for 
7 = 2 may be identical to the one of a pressureless fluid in GR (e.g. CDM). 
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6 Evolution of sub-horizon perturbations 



The results of our background tests do not impose negative values for a with a high confidence 
level. This could mean that the GCG has indeed a = meaning that the dark sector of the 
universe is actually an ordinary adiabatic fluid with a small, constant and negative pressure. 
Hence, we can argue that at high redshifts the GCG does not look exactly like standard 
CDM {p,dm = 0). 

Our goal in this section is to revisit the issue concerning the equivalence of GCG with 
a = and the ACDM model. Obviously both cases produces the same background expansion. 
However, does this equivalence holds at first order in the perturbations? 

We present an analysis of the Meszaros effect which describes the formation of dark 
matter structures during the first stages of the matter dominated epoch. We assume a 
single fiuid description, i.e. we neglect the contribution of baryons and radiation. This is a 
reasonable approximation if we want to track the growth of sub-horizon perturbations during 
the matter dominated epoch. For the technical details, we follow the analysis of [47]. 

For the linear perturbations of (5.18) we define the velocity scalar which is associated 
with the peculiar velocity by Su^^ = —kv/a. At first order, the (0-0) component of the 
Einstein equation reads 

- k"^^ - m^' - = ^n^s , (e.i) 

where 6 is the density contrast. During the time interval between the kinetic decoupling of 
DM particles from the primordial plasma and the epoch of matter radiation equality {zeq ~ 
3300 for the ACDM model) the sub-horizon DM perturbations grow only logarithmically 
with the scale factor. After Zeq the DM perturbations obey to 5 oc a. This is the main result 
behind the so called Meszaros effect [49]. 

We study here an example in which the GCG with a = behaves differently from the 
ACDM. We will show that the growth of sub-horizon GCG matter perturbations could be 
different from the standard CDM. In order to obtain a Meszaros-like equation for the GCG 
we make use of the covariant conservation of the energy-momentum tensor (T^^ = 0). At 
first-order we find 

S' = -m5 [cl - Wc) - (1 + Wc) {kv - W) = , (6.2) 

and 

v' = -n{l- 3-u;,) V + k^ + —^5 = , (6.3) 

for the energy and momentum balances, respectively. 

For small scales we can neglect the term ^' in Eq. (6.2) and we also take the sub-horizon 
limit of the Poisson equation Eq. (6.1). Together with Eq. (6.3) we find a Meszaros-like 
equation for the GCG with a = 0, i.e. = 0: 



da'' + 



3 15 2 'iawc dH dwr. 
—Wc + 9w,. — ; 3a- 



2 2 H da da 



dSqcq 

a — ^ + 
da 

5gcg = 0. (6.4) 
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The standard equation for CDM 



a' 




+ 3 + 



a dH] ddcdm 3 



(6.5) 



da? 



H da da 2 



can be obtained from (6.4) if we set Wc — 0, resulting in the expected solution Scdm ^ ^ 
during the matter dominated phase. 

Obviously we note that GCG sub-horizon matter perturbations do not evolve in the 
same way as standard pressureless matter present in the ACDM model. Thus, we provide an 
example where the models have different predictions contradicting Ref. [32]. 

We show in Figure 9 the evolution of a typical sub-horizon scale after the equality 
{zeq ~ 3000). We assume that both CDM and GCG have the same initial conditions at that 
time. We set them by generating the power spectrum at that time for the scale k = 0.2 h 
Mpc~^ with help of the CAMB code [51]. We then solve numerically Eq. 6.4 for A = 0.05, 0.4 
and 0.90. The perturbations in the GCG fluid are strongly suppressed in comparison to the 
CDM case. They do not reach the the non-linear regime 5 ~ 1. This means that a small 
negative pressure of the matter fluid at early times drastically prevents structure formation. 
The small scales would be the most affected ones and only numerical simulations would 
provide the correct matter distribution. However, semi-analytical studies touching the non- 
linear regime [10, 48] could supply additional information about the final clustering stage. 

7 Conclusions 

We performed a Bayesian analysis of the background behaviour of the GCG model using 
H[z), SNIA, CMB and BAO datasets. We focused particularly on the parameter a, on 
which a huge literature already exists. Our result is that a = — O.OSTtoJas, at the 2a level, 
i.e. negative values of a seem to be favoured over the positive ones. The uncertainty is too 
large for us to claim that a "tension" with perturbative tests (which constrain |a| < 10~^ 
at 2(t) does exist, but we speculate on the consequences of this occurrence, which is not far 
from being settled, given the ongoing and forthcoming observational programs which collect 
more precise data day after day. That kind of inconsistency would immediately rule out the 
GCG model, since it is unacceptable to obtain two different set of parameters depending on 
the observational test applied. We assume this to occur and figure out how to possibly save 
the GCG unification paradigm. We introduce a fiuid model in Rastall's theory which exactly 
behaves as the GCG at the background level (i.e. it provides the same evolution for H{z)) 
and discuss its perturbative properties. It turns out that small and positive speeds of sound 
are compatible with a < 0. On the other hand, we show that the evolution of perturbations 
is ruled by an effective speed of sound which is identical to the GCG one in GR, as if Rastall's 
theory introduced a sort of "geometric entropy". Therefore, as one should expect, the results 
are catastrophic, with the density contrast and the gravitational potential growing too fast 
(for a < 0) for being in agreement with observation. As a possible salvation, we indicate a 
scalar field description of the GCG in Rastall's theory, where the background expansion and 
the effective speed of sound can be fixed independently (and the latter to zero). 

We also address the issue of the a — )■ limit, which is considered indistinguishable from 
the ACDM model. We show that this is not the case. The small negative pressure of the 
GCG affects the evolution of the GCG density contrast during the matter dominated phase, 
where it is expected to behave as CDM and provides a different evolution for the sub-horizon 
matter perturbations. Indeed, this analysis challenges the matter epoch of the GCG showing 
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that sub-horizon structures cannot form as in the standard ACDM case. Thus, it could be 
very interesting to perform numerical simulations for matter fluids with a hydrodynamical 
evolution distinct from the standard pressureless case. In particular, a fluid with a constant 
negative pressure. This analysis will surely provide a different clustering pattern at galactic 
scales. 
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Figure 1. One-dimensional PDF for the parameter a under the restriction a > —1 using H{z) (upper 
left panel), CMB (upper right panel), SNIa (lower left panel) and BAO (lower right panel). We have 
used here delta priors of h. 
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Figure 2. One-dimensional PDF for the parameter A under the restriction a > — 1 using H{z) 
(upper left panel), CMB (upper right panel), SNIa (lower left panel) and BAO (lower right panel). 
We have used here delta priors of h. 



Figure 3. One-dimensional PDF for the parameter a under the restriction a > — 1 using li{z) (upper 
left panel), CMB (upper central panel), SNIa (upper right panel), BAO (lower left panel) and the 
combination of these four tests (lower right panel). We have used a flat prior for h and marginalized 
over it. 
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Figure 4. One-dimensional PDF for tlie parameter A under the restriction a > — 1 using H{z) 
(upper left panel), CMB (upper central panel), SNIa (upper right panel), BAO (lower left panel) and 
the combination of these four tests (lower central panel). In the lower right panel, we present the 
contour plots (1, 2, 3 a) for the total PDF in the (a, A) plane. We have used a flat prior for h and 
marginalized over it. 



Figure 5. One-dimensional PDF for the parameter a, with no restriction, using H{z) (upper left 
panel), CMB (upper central panel), SNIa (upper right panel), BAO (lower left panel) and the combi- 
nation of these four tests (lower right panel). We have used a flat prior for h and marginalized over 
it. 
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Figure 6. One-dimensional PDF for tlie parameter A, with no restriction on a, using H{z) (upper 
left panel), CMB (upper central panel), SNIa (upper right panel), BAO (lower left panel) and the 
combination of these four tests (lower central panel). In the lower right panel, we present the contour 
plots (1, 2, 3 cr) for the total PDF in the (a. A). We have used a flat prior for h and marginalized 
over it. 
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Figure 7. Evolution of the equation of state parameter w and of the speed of sound given in 
Eqs. (5.11) and (5.12) as functions of the scale factor. The parameter 7 has been chosen 7 = 0.8 
(black solid lines), 7 = 1 (i.e. the GR limit, red dashed lines) and 7 = 1.1 (blue dash-dotted line). 
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Figure 8. Evolution of $ and 5 as functions of the scale factor. The GCG parameters have been 
fixed as a = —0.1 and A = 0.7. The parameter 7 has been chosen 7 = 0.5 (black solid lines), 7 = 1 
(i.e. the GR limit, red dashed lines) and 7 = 2 (blue dash-dotted line). For comparison, the ACDM 
lines are depicted as the magenta dotted ones. 
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Figure 9. Evolution of the sub-horizon density contrast 5. The standard prcssurclcss CDM {p = 0) 
is shown as the solid line. For the GCG with a = we show the evolution of 5 for three different 
values of A, as indicated in the figure. 



